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A new jet observable, dipolarity, is introduced that can distinguish whether a pair of subjets arises 
from a color singlet source. This observable is incorporated into the HEPTopTagger and is shown to 
improve discrimination between top jets and QCD jets for moderate to high pr- 



Introduction 

The impressive resolution of the ATLAS and CMS de- 
tectors means that a typical QCD jet at the LHC deposits 
energy in 0(10—100) calorimeter cells. Such fine-grained 
calorimetry allows for jets to be studied in much greater 
detail than previously, with sophisticated versions of cur- 
rent techniques making it possible to measure more than 
just the bulk properties of jets {e.g. event jet multiplic- 
ities or jet masses). One goal of the LHC is to employ 
these techniques to extend the amount of information 
available from each jet, allowing for a broader probe of 
the properties of QCD. The past several years have seen 
significant progress in developing such jet substructure 
techniques. A number of general purpose tools have been 
developed, including: (i) top-tagging algorithms designed 
for use at both lower [1, 2 and higher [3^ pt as well as 
(ii) jet grooming techniques such as filtering [4], pruning 
[5 , and trimming [6 , which are designed to improve jet 
mass resolution. Jet substructure techniques have also 
been studied in the context of specific particle searches, 
where they have been shown to substantially extend the 
reach of traditional search techniques in a wide variety of 
scenarios, including for example boosted Higgses [HIT], 
neutral spin-one resonances [8], searches for supersym- 
metry [9j, and many others p!QHT4] . Despite these many 
successes, however, there is every reason to expect that 
there remains room for refinement of jet substructure 
techniques. 

Top tagging algorithms have reached a mature level of 
development in recent years [Ir^ 15 - 20^ . A variety of dif- 
ferent algorithms employ primarily kinematic observables 
like jet masses and the helicity angle in increasingly 
sophisticated ways, allowing for efficient discrimination 
between top jets and ordinary QCD jets. One direction 
that has received less attention is the use of observables 
that do not map onto the kinematics of hard partons. 
This article introduces a new non-kinematic observable 
that can be used in top tagging to gain additional back- 
ground rejection and should have applications outside of 
top tagging. 

A distinguishing aspect of hadronic top decays is that 



the jets from the decay belong to a color singlet 
configuration. For prw ^ these jets become close 
together and will often be clustered within a single jet. 
The radiation pattern of the decay products, which 
is controlled by the color configuration, has a distinctive 
form with most of the radiation clustered in between the 
two jets. This QCD analog of the Chudakov effect of- 
fers an additional handle for top discrimination on top of 
kinematic observables. 

The organization of this article is as follows. First the 
concept of color flow is briefly reviewed, and connection 
is made to the jet observable pull [21 . Next the pro- 
posed "dipolarity" observable is introduced and explored 
in a general context. Dipolarity is then incorporated into 
the HEPTopTagger, and the performance of the modified 
tagger is tested on Monte Carlo event samples. Finally, 
there is some discussion of the results as well as possible 
applications of dipolarity beyond top tagging. 
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FIG. 1: Legoplot for a top jet reconstructed by the 
HEPTopTagger, with hard substructure identified by a com- 
bination of filtering and a fractional mass-drop criterion. The 
orange cells correspond to soft radiation associated with the 



Color flow and pull 

Within the context of top tagging, several jet observ- 
ables have been defined that go beyond the kinematics 
of hard partons. These include a number of jet shape 
observables such as spherocity [18] , planar fiow [20l |22] , 
A^-subjettiness [23], and template overlap [24 . The jet 
observable defined in the next section draws from the 
complimentary information offered by color fiow. In a 
QCD event, radiation is controlled by the kinematics of 
the hard partons as well as by how color indices are con- 
tracted together (color fiow) . Partons whose color indices 
are contracted together are color-connected, with a color 
string stretching between the two color sources. For ex- 
ample, the two quarks in the hadronic decay of a color 
singlet like the Higgs form a color dipole whose radia- 
tion pattern is contained primarily within a pair of cones 
around the two quarks, with a tendency for more radia- 
tion to occur in the region between the two quarks [25 . 

Color fiow arguments of this sort have motivated at- 
tempts to use QCD radiation patterns for event dis- 
crimination, e.g. mini-jet vetoes in Higgs searches [26^. 
More recently, the authors of [21] introduced a jet ob- 
servable dubbed pull, which is a pT-weighted vector in 
rapidity-phi space that is constructed so as to point from 
a given jet to its color-connected partner (s). Although 
pull has been shown to offer some discrimination in par- 
ticle searches [27 , it does not seem well-suited to tagging 
boosted hadronic tops. The most straightforward way to 
incorporate pull into a top tagging algorithm is to mea- 
sure the pull of two subjets that reconstruct the and 
check whether each subjet's pull vector points towards 
the other subjet. A problem with this approach is that 
the pull vectors are sensitive to how the jet is broken 
down into two subjets. For a lopsided distribution of the 
into two subjets, one of the subjets will consist of 
only a small handful of calorimeter cells, and as a conse- 
quence its pull will be sensitive to statistical fiuctuations 
and contamination. Even for a broken down into 
two subjets more symmetrically, the pull vectors can de- 
pend sensitively on the precise boundary drawn between 
the two subjets, which itself is a noisy function of the 
particular jet clustering algorithm being used. Antici- 
pating the potential difiiculties of incorporating pull into 
a top tagging algorithm, we explore an alternative ap- 
proach in which the entire radiation pattern of the 
is considered simultaneously. This simple idea leads us 
to jet dipolarity, which we now define. 

Dipolarity 

Consider a jet, J, with two subjets, ji and j2, whose 
centers are located at pseudorapidities rji and 7^2 and 
azimuthal angles 0i and ^2, respectively. For each 
calorimeter cell (r^^, (pi) with transverse momentum pTi 
let Ri be the (minimum) euclidean distance in the rj-cj) 
plane between (t^^, and the line segment that runs 



from (771, ^1) to (7^2, ^2)- Dipolarity is defined as the 
PT-weighted sum 

where = (^1 — ^2)^ + (0i — 02)^- Dipolarity is an 
infrared and collinear (IRC) safe observable so long as 
the algorithm used to identify J, ji and j2 is IRC safe. 
Notice that dipolarity, which is essentially a two-subjet 
observable, requires the centers of ji and ^2 as input, 
although it does not require that the constituents of J be 
partitioned between ji and ^2. The centers of ji and j2 
can be determined by whatever procedure is convenient 
for the particular application. For example one could 
choose the centers of ji and so as to minimize the sum 
in 

Dipolarity will be small when most of the radiation 
within the jet J occurs in the region between the two 
subjets ji and ^2 and will be large whenever a substantial 
amount of radiation is found elsewhere. As a consequence 
of the weighting with respect to in ([T]), V receives 
large contributions from semisoft radiation away from the 
cores of ji and ^2. It is this semisoft radiation away 
from the cores of ji and j2 that is expected to refiect the 
color configuration of J. The weighting in ([T]) does not 
know about the exact radiation pattern of a color singlet; 
nevertheless, we expect that color singlets that decay into 
two jets will have small P, while radiation emitted by 
colored objects will tend to yield larger values of V. 

This expectation can be fieshed out more explicitly by 
considering the emission pattern of a third parton with 
energy uj from a pair of partons in a particular color 
configuration, see e.g. [25]. In the eikonal approximation 
(cj 0) one finds the radiation function for a color sin- 
glet to be 

WM^)^ , ^ X'^, ^ ^. (2) 

while for two partons color-connected to the beam we 
have instead 

dr]d(j) 

'1 ^beam ) 

where 

X(7/, ^; r]i, (t)i) = cosh{r] - rji) - cos((/) - (/)i) (4) 

The resulting radiation patterns are depicted in FIGj2j 
One sees explicitly that the color singlet has its radiation 
clustered in the region between the two partons, whereas 
for partons color-connected to the beam, a substantial 
amount of radiation is emitted towards the beam. Using 
the expressions in ([2| and ([3| to calculate V gives the pre- 
diction Vns 2Ps; although this is approximately what 
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FIG. 3: Schematic for a collection of QCD jets whose kine- 
matics fake the top. The upper figures show various possi- 
bilities for quarks and gluons that undergo two branchings. 
The bottom figures show the corresponding large Nc color di- 
agrams, with dipole radiation patterns superimposed across 
color dipoles. Only the rightmost color configuration, which 
is suppressed by factors of Ca/Cf with respect to the others, 
matches the radiation pattern of an actual top. 

ent color configurations fakes the kinematics equally well. 
The dipolarities of the subjets, however, will be broadly 
distributed in accord with their different color configu- 
rations. For instance, gluon jets are known to give the 
largest fake rates for top jets as a consequence of their 
larger Casimirs which more often result in wide angle 
branchings with significant mass drops. FIG. [3] illustrates 
how gluon jets, with their distinct color configurations, 
radiate differently from top jets. All of this suggests that 
the dipolarity of the in a hadronic top decay is well- 
suited as a discriminant in top tagging algorithms. 



FIG. 2: Top: Eikonal radiation pattern dpr/dridcj) for a color 
singlet with Ai?=0.9, typical for a originating from a 
top with pt ^ 300 GeV. Bottom: As above with the partons 
instead color-connected to the beam (left/right-going parton 
connected to the left/right beam). Contours differ by powers 
of e. For the color singlet the radiation is mostly found in 
the region between the two subjets. For the background- like 
color configuration, the radiation is pulled towards the beam. 
Note that an absolute comparison cannot be made between 
the figures, since the collinear singularities in ^ and ([3| are 
not regulated. 



is found from Monte Carlo calculations, expressions ([2| 
and (|3| do not yield dipolarity distributions in quantita- 
tive agreement with the Monte Carlo. Given the crude- 
ness of the approximations that went into these expres- 
sions, this discrepancy is not surprising; a more accurate 
estimate of V for various color configurations could be 
obtained by using antenna patterns as in [28 . 

Dipolarity can be used within the context of top tag- 
ging to reduce QCD backgrounds. Consider a collection 
of fat QCD jets originating from parton branchings with 
identical kinematics but different color configurations as 
illustrated in FIG. [S) If one of the QCD jets fakes the 
kinematics of a top quark decay, then each of the differ- 



HEPTopTagger 

To test whether dipolarity makes an effective dis- 
criminant, cuts on dipolarity are incorporated into the 
HEPTopTagger (TJ |2], which is designed to work effec- 
tively at intermediate boost, with 200 GeV ^ Pt ^ 
800 GeV. The high efficiency of the HEPTopTagger at 
these Pt makes it a good candidate for such a modifi- 
cation because dipolarity cuts are expected to be most 
effective at intermediate pt- This is because at lower pt 
contamination from pile-up and the underlying event be- 
comes more of a concern as the top jets become fatter 
and fatter, while at higher pt the finite resolution of the 
detector makes it difficult to get an accurate handle on 
radiation patterns. Furthermore, the multibody filtering 
implemented by the HEPTopTagger results in accurate re- 
construction of the W^. The HEPTopTagger algorithm 
is defined as follows: 

1. Using the Cambridge /Aachen algorithm cluster the 
event into fat R = 1.5 jets. 

2. Break each fat jet j into hard subjets using the 
following mass-drop criterion: undo the last stage 
of clustering to yield two subjets ji and j2 (with 
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rrij^ > rrij^), keeping both ji and j2 if rnj^ < O.Srrij 
and otherwise dropping j2; repeat this procedure 
recursively, stopping when the rrij. drop below 30 
GeV. 

3. Consider in turn all possible triplets of hard sub- 
jets. First, filter^ each triplet with a resolution 
^filter = niin(0.3, Ai?^j/2). Next, using the five 
hardest constituent subjets of the filtered triplet 
calculate the jet mass mfiit- Finally, choose the 
triplet whose mmt lies closest to m^. 

4. Recluster the five filtered constituents chosen in 
step 3 into exactly three subjets ji, j2, and js or- 
dered in descending pt- Accept the fat jet as a 
top candidate if it passes any of the following three 
pairs of mass cuts: 
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5. Finally, require that the total pr of the three sub- 
jets defined in step 4 be greater than 200 GeV.^ 

Dipolarity cuts are introduced into the HEPTopTagger 
by modifying step 4 above. For a top candidate that has 
passed one of the three pairs of mass cuts we calculate 
the dipolarity of the as identified by the mass cut: 
e.g. for a top candidate that satisfies ii) and ii') the is 
identified as ji + j2- If more than one of the pairs of mass 
conditions is satisfied in step 4, we choose the smaller 
dipolarity. We find that this procedure performs better 
than calculating the dipolarity of the pair of subjets that 
reconstructs mw most accurately. 

In addition to introducing dipolarity cuts, we also 
make cuts on the filtered mass of the reconstructed top, 
mfiit, which is not done in the original HEPTopTagger, 



^ That is, recluster the constituents of the triplet using the Cam- 
bridge/Aachen algorithm with jet radius i^fiiter- Filtering was 
introduced in 4 . 

The HEPTopTagger does not make use of 6-tagging, which is a 
natural extension to the algorithm that can result in significant 
improvements in background rejection. Since dipolarity cuts are 
orthogonal to b-tagging, we do not explore the use of b-tagging 
in this paper. 
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FIG. 4: Dipolarity distributions for W^s reconstructed by 
the HEPTopTagger and passing default mass cuts with mfiit G 
[150 GeV, 210 GeV]. Thick solid lines indicate central values, 
whereas thin dashed lines correspond to values at 10% and 
90%. Here and throughout the pr is that of the fat R — 1.5 
jet. For all pr the central value of the dipolarity for the back- 
ground is (9(50% — 100%) larger than for the signal. This fig- 
ure uses the HERWIG event samples; the PYTHIA event samples 
yield similar distributions. The dot-dash line roughly indi- 
cates where dipolarity cuts are made at the S=20% working 
point. 



where the cuts have been chosen so as to avoid any ex- 
plicit mass scales. We introduce cuts on mfiit for two 
reasons. The first is to improve background rejection. 
The second and main reason is that we are interested 
in determining whether dipolarity cuts are essentially or- 
thogonal to cuts on kinematic observables. To do this we 
must ensure that the HEPTopTagger is using a full com- 
pliment of kinematic cuts, including cuts on mfiit- In a 
particular application, cuts on mfiit may be undesirable. 
In that case, the inclusion of dipolarity cuts would result 
in a larger improvement of background rejection. 

Note that the ji selected in step 4 contain only the hard 
substructure of the fat jet. Some amount of soft radiation 
has been thrown out by filtering and the mass drop crite- 
rion. To eflFectively gauge whether the radiation pattern 
of the reconstructed is consistent with the expected 
dipole radiation pattern, it is important to include some 
of the discarded soft radiation. We find that the crite- 
rion used to select the radiation included in calculating 
the dipolarity of the has significant impact on the 
ultimate utility of dipolarity as a discriminant. In par- 
ticular diflFerent criteria will lead to dipolarity distribu- 
tions that are more or less correlated with the kinematic 
observables considered by the HEPTopTagger. Applying 
dipolarity in another context would likely require this 
criterion to be carefully reworked so as to maximize per- 
formance. In the present case we find that the following 
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criterion, which aims to capture as much of the radiation 
emitted by the as possible, while minimizing possible 
contamination, to be most effective. In addition to the 
hard radiation from the two subjets, we include all 
soft radiation contained within the pair of cones centered 
around the two hard subjets, fixing the radius of the 
cones to be Ai?/V^, where AR is defined between the 
two hard subjets. Furthermore we exclude any ra- 
diation contained within the smallest cone that encloses 
the hard b subjet. Note that angular ordering implies 
that the majority of the radiation emitted by the 
is within the pair of cones of radius AR. We choose our 
cones to be somewhat smaller to minimize contamination 
from underlying event /pile-up as well as the b subjet. See 
FIG. [T]to see this selection criterion at work on a sam- 
ple top jet. The orange cells in the figure correspond to 
radiation that has been discarded by mass-drop filtering 
but that is included in calculating the dipolarity of the 
candidate W^. 



Testing Dipolarity 

In order to make a meaningful comparison between 
the performance of the HEPTopTagger with and without 
dipolarity cuts, it is not enough to leave the kinematic 
cuts employed by the HEPTopTagger fixed at their default 
values. Instead it is important to optimize the cuts to 
yield the largest background rejection at each given signal 
efficiency. This optimization is performed by a custom 
Monte Carlo code that finely samples the space of cuts. 
Specifically, a scan is performed at discrete values of i^min 

70% mw/'mt < Rmin < 98% mw/'mt 



with step sizes of 1% mw/'mt' For each value of 

i^min, ^max IS choSCU tO be i^max 2.0 X TTlw / rUt - 

Rmin- In addition we simultaneously optimize over cuts 
"mit jxiin < ^fiit < max, ^ < ^max, and, in step 
4 of the HEPTopTagger, /cut < arctan mis/mi2 and 
^cut ^ ^23/^123- The remaining cuts in step 4 of the 
HEPTopTagger are left at their default values, since these 
are less important for background rejection. Addition- 
ally, the two parameters that define the mass drop cri- 
terion remain fixed at their default values. The hard 
substructure cutoff of 30 GeV is sensitive to detector ef- 
fects, which can only be crudely mocked-up without a 
full detector simulation. See TABLE [l| for a sample set 
of cuts. 

We use three different event samples for evaluating the 
performance of the modified HEPTopTagger. These event 
samples (with center of mass energy of 7 TeV) belong to 
a set of benchmark event samples that have been made 
publicly available by participants of BOOST 2010 [29]. 
The first event sample is generated by HERWIG 6.510 [30^ 
with the underlying event simulated by JIMMY [31 , which 
has been configured with a tune used by ATLAS. The 
second is generated by PYTHIA 6.4 [32^ with Q^-ordering 



and the 'DW tune for the underlying event. The third 
is generated by PYTHIA 6.4 |32J with pT-ordering and the 
'Perugia' tune for the underlying event. See [33] for more 
details. For signal jets we use the hardest jet in each event 
of a Standard Model hadronic tt sample, excluding jets 
with 1 7^ I > 2.5. For background jets we use the hardest 
jet in each event of a Standard Model dijet sample, again 
excluding jets with |7^| > 2.5. For jet clustering we use 
the Cambridge- Aachen (CA) algorithm [34l |35| as im- 
plemented by Fast Jet 2.4.2 [36 . In order to simulate 
the finite resolution of the ATLAS or CMS calorimeters, 
particles in each event are clustered into 0.1 x 0.1 cells in 
r] — (j) space and then combined into massless four- vector 
pseudoparticles that are fed into Fast Jet. We have also 
checked that imposing a low energy cutoff of 1 GeV on 
each cell results in only a mild degradation of background 
rejection. 
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FIG. 5: Background mistag rates with and without dipolarity 
cuts at a fixed signal efficiency of 20% for HERWIG (top) and 
PYTHIA (bottom) . The mistag rate at a given pro is calculated 
from a pr window of 100 GeV centered at pro- Note that, 
as a consequence, each point is not statistically independent. 
Error bands are statistical. 

FIG. [4] shows the dipolarity distributions that result 
from the modified HEPTopTagger for a default set of 
mass cuts. As expected the dipolarity of the recon- 
structed in top jets is smaller than for the QCD 
background. FIG. [5] shows the improvement in perfor- 
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Operating Point 


Icut 


^cut 


^min 


^max 


min 


'^t max 


B{%) 


Low pt without V cut 


0.45 


0.41 


0.92 




159 GeV 


195 GeV 


0.47 


0.48 


0.41 


Low Pt with V cut 


0.37 


0.39 


0.80 


0.021 


154 GeV 


199 GeV 


0.41 


0.38 


0.30 


High Pt without V cut 


0.47 


0.40 


0.93 




158 GeV 


199 GeV 


0.92 


0.79 


0.59 


High Pt with V cut 


0.36 


0.38 


0.88 


0.023 


154 GeV 


196 GeV 


0.58 


0.58 


0.39 



TABLE L Sample optimized operating points at 5* = 20% based on an equal admixture of all three event samples for maximum 
statistics. The resulting background mistag rates (B) are shown for each of the three event samples with HERWIG, PYTHIA 'DW, 
and PYTHIA 'Perugia' arranged from left to right. Including dipolarity cuts loosens mass cuts while improving background 
rejection. The low pt samples have 200 GeV <Pt< 500 GeV while the high pt samples have 400 GeV <Pt< 800 GeV. 




Signal efficiency (%) 

FIG. 6: Signal efficiency vs. background mistag rate for the 
HEPTopTagger with pt G [400 GeV, 500 GeV] and HERWIG 
event samples. At lower signal efficiencies the inclusion of 
dipolarity cuts leads to a sizeable improvement in background 
rejection. Statistical error bars, which are a relative 25% at 
the lowest mistag rate, are not shown. 



mance that results from including dipolarity cuts into 
the HEPTopTagger as a function of pt at a fixed signal 
efficiency of 5* = 20%. Finally, FIG. [6] shows how, for 
Pt ^ [400 GeV, 500 GeV], dipolarity cuts improve back- 
ground rejection at signal efficiencies S ^ 35%. 

At a fixed signal efficiency oi S = 20%, dipolarity 
cuts lead to a sizable decrease in the mistag rate for 
350 GeV <pt ^800 GeV with the largest decrease at in- 
termediate For the HERWIG event samples the mistag 
rate decreases by as much as ~ 50%, whereas for the 
PYTHIA samples the mistag rate decreases by a more 
modest amount, ~ 30%. Differences in the underlying 
event (UE) model do not explain this disagreement; for 
instance, repeating the analysis with the UE turned off 
results in background mistag rates that are only some- 
what smaller, with both HERWIG and PYTHIA dropping 
by similar amounts. This suggests that the difference 
between the HERWIG and PYTHIA mistag rates most likely 
arises from the parton shower. 



Discussion 

This article has introduced a new jet substructure ob- 
servable that is useful in discriminating among different 
color configurations in jets that have large mass drops. 
This discrimination is of interest, since such jets often 
arise from decays of boosted heavy particles. Incorpo- 
rating this discriminant into a top tagging algorithm re- 
sults in QCD background mistag rates that are lower 
by as much as 50%; the exact mistag rate, however, 
shows considerable sensitivity to the details of the par- 
ton shower. Specifically HERWIG event samples result in a 
larger improvement in background rejection than is found 
for PYTHIA. We suspect that HERWIG, which uses angu- 
lar ordering, does a better job of simulating the effects 
of color coherence than PYTHIA, which uses or pt- 
ordering in combination with angular vetoes. This could 
explain why the dipolarity of the is a more discrim- 
inating observable in the case of the HERWIG event sam- 
ples. With an understanding of the origin of this dif- 
ference, comparisons to measurements at the LHC could 
help improve the description of QCD radiation. It would 
be interesting to understand how this difference arises 
from the details of the parton shower; doing so, however, 
lies outside the scope of this paper. 

Validating how well color flow effects as modeled by 
Monte Carlo event generators match what is observed 
in collider experiments is only beginning to be studied 
actively. Understanding color flow in detail is a diffi- 
cult problem; for example, QCD predictions for radia- 
tion patterns can be affected by non-global logarithms, 
see e.g. [40 . Therefore validating theoretical predictions 
against data will be critical in reducing the theoreti- 
cal uncertainty associated with how dipolarity and other 
color flow observables are modeled by Monte Carlo calcu- 
lations. A few color coherence studies performed at the 
Tevatron showed spatial correlations between the third 
and second hardest jets in pp collisions, and HERWIG was 
shown to provide a better description of the data than 
PYTHIA More recently, the color of the in ti 

events was studied, and agreement between theory pre- 
dictions for jet pull and data was shown [38] . 

Jet dipolarity should be useful in a broader set of ap- 
plications to both Standard Model and beyond the Stan- 
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dard Model physics. Possible directions for future re- 
search include: (i) applications of dipolarity to a collider 
search for heavy color singlet resonances that decay to 
tt; (ii) applications to standard model W^/Z^ physics; 
(iii) applications to heavy color singlet resonances that 
decay to or Z^Z^; (iv) applications to cascade 

decays of super symmetric particles; (v) inclusion of dipo- 
larity into other top-tagging algorithms; (vi) applications 
to the decay of new particles into novel color configura- 
tions such as in the decay of the LSP in supersymmetric 
models with baryonic R-parity violation; and (vii) modi- 
fying V to more closely correspond to the exact radiation 
pattern of a color singlet. Each of these directions could 
make an interesting laboratory for further development 
of jet substructure techniques. 
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